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Abstract 

An instance of the maximum mixed graph orientation problem consists of a mixed graph and a col- 
lection of source-target vertex pairs. The objective is to orient the undirected edges of the graph so as to 
maximize the number of pairs that admit a directed source-target path. This problem has recently arisen 
in the study of biological networks, and it also has applications in communication networks. 

In this paper, we identify an interesting local-to-global orientation property. This property enables us to 
modify the best known algorithms for maximum mixed graph orientation and some of its special structured 
instances, due to Elberfeld et al. (CPM '11), and obtain improved approximation ratios. We further proceed 
by developing an algorithm that achieves an even better approximation guarantee for the general setting of 
the problem. Finally, we study several well-motivated variants of this orientation problem. 
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1 Introduction 



An instance of the maximum mixed graph orientation problem consists of a mixed graph G = {V, Ed U E\j) 
with n vertices, such that Ed and E-\j indicate the sets of directed and undirected edges, respectively. An 
additional ingredient of the input is a collection P C F x F of source-target vertex pairs. A source-target 
vertex pair {s, t) £ P is called a request. The objective is to orient G in a way that maximizes the number of 
satisfied requests. An orientation of G is a directed graph G = (V, Ed U Etj), where E\j is a set of directed 
edges obtained by choosing a single direction for each undirected edge in E\]. A request (s, t) is said to be 
satisfied under an orientation G if there is a directed path from s to t in G. 

One may assume without loss of generality that the mixed graph G is acyclic, that is, a graph that has 
no cycles. This assumption holds since any instance of maximum mixed graph orientation can be reduced to 
another instance in which the underlying mixed graph is acyclic without affecting the number of requests that 
can be satisfied [23, 10]. Indeed, if the input graph contains cycles, one can sequentially contract them one 
after the other. In each step, the undirected edges of an arbitrary cycle are all oriented in the same direction. 
In particular, if this cycle contains directed edges then the undirected edges are oriented in a consistent way 
with those edges. As a result, every pair of vertices on this cycle admits a directed path between them, and 
thus, the cycle can be contracted. One can easily validate that the resulting mixed acyclic graph consists of 
undirected components, each of which must be an undirected tree, and those components are connected by 
directed edges in a way that does not produce cycles. The maximum mixed graph orientation problem draws 
its interest from applications in network biology and communication networks: 

Network biology. Recent technological advances, such as yeast two-hybrid assays [12] and protein co- 
immunoprecipitation screens [16], enable detecting physical interactions in the cell, leading to protein-protein 
interaction (PPI) networks. One major caveat of those PPI measurements is that they do not reveal informa- 
tion about the directionality of the interactions, namely, the directions in which the signal flows. Since PPI 
networks serve as the skeletons of signal transduction in the cell, inferring the hidden directionality informa- 
tion may provide insights to the inner working of the cell. Such an information may be inferred from causal 
relations in those networks [24]. One such source of causal relations is perturbation experiments, in which 
a gene is perturbed (cause) and as a result, other genes change their expression levels (effects). A change of 
expression of a gene suggests that the corresponding proteins admit a path in the network, and in particular, it 
is assumed that there must be a directed path from the causal gene to the affected gene. 

Up until this point in time, the above-mentioned scenario can be modeled as a special instance of the 
maximum mixed graph orientation problem in which one is interested to orient the edges of an undirected 
network in a way that maximizes the number of cause-effect pairs that admit a directed path from the causal 
gene to the affected gene. However, in the more accurate biological variant, there are several interactions 
whose directionality is known in advance. For instance, protein-DNA interactions are naturally directed from 
a transcription factor to its regulated genes, and some PPIs, like kinase-substrate interactions, are known to 
transmit signals in a directional fashion. Therefore, in general, the input network is a mixed graph. 

Communication networks. A unidirectional communication network consists of communication links that 
allow data to travel only in one direction. One main benefit of such communication links is that the data of 
the device on one side is kept confidential while it may still access the data of the device on the other side. 
As a consequence, unidirectional networks are most commonly found in high security environments, where a 
connection may be made between devices with differing security classifications. For example, unidirectional 
communication links can be used to facilitate access to a vulnerable domain such as the Internet to devices 
storing sensitive data. The maximum mixed graph orientation problem captures the interesting scenario in 
which one is interested to design a unidirectional network that maximizes the number of connection requests 
that can be satisfied in a secure way. We remark that unidirectional networks have also been studied in 
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distributed and wireless ad hoc settings (see, e.g., [2, 1, 21] and the references therein), where a common 
focus is on algorithmic questions that arise in a given unidirectional network. Here, we are rather interested 
in the question of how to design such a network while optimizing some performance guarantees. 

1.1 Previous work 

Arkin and Hassin [3] seem to have been the first to study the problem of orienting mixed graphs. They focused 
on the decision problem con^esponding to maximum mixed graph orientation, and demonstrated that it is NP- 
complete. Elberfeld et al. [10] observed that the reduction in their proof implies that the maximum mixed 
graph orientation problem is NP-hard to approximate to within a factor of 7/8. Silverbush, Elberfeld, and 
Sharan [23] devised a polynomial-size integer linear program formulation for this problem, and evaluated its 
performance experimentally. Recently, Elberfeld et al. [10] developed several polylogarithmic approximation 
algorithms for special instances of the problem in which the underlying graph is tree-like, e.g., when the graph 
has bounded treewidth. In addition, they developed a greedy algorithm for the general setting that achieves 
Q,{\/{M'^ log n))-approximation, where M = max{n, \P\} and c = l/\/2 0.7071. 

Medvedovsky et al. [22] initiated the study of the special setting of maximum graph orientation in which 
the underlying graph is undirected, that is, when there are no pre-directed edges. They proved that it is 
NP-hai^d to approximate this problem to within a factor of 12/13, even when the graph is a star. They also 
proposed an exact dynamic -programming algorithm for the special case of path graphs, and a 0(1/ log n)- 
approximation algorithm for the general problem. Gamzu, Segev and Sharan [15] utilized the framework 
developed in [14] to obtain an improved 0(log log n/ log n) -approximation ratio (see also [9]). Very recently, 
Dorn et al. [8] studied this problem from a parameterized complexity point of view. They presented several 
fixed-parameter tractability results. Further research focused on other variants of this undirected orientation 
problem. For example, Hakimi, Schmeichel, and Young [17] studied the special setting in which the set of 
requests contains all vertex pairs, and developed an exact polynomial-time algorithm. 

1.2 Our results 

We identify a useful structural property of requests crossing through a junction vertex. Informally, this prop- 
erty guarantees that if a set of requests is locally satisfiable then it can also be satisfied globally. Using this 
property, we can slightly modify the algorithms developed by Elberfeld et al. [10], and obtain improved ap- 
proximation ratios. For example, we eliminate a logarithmic factor from their polylogarithmic approximation 
ratio for the case that the underlying graph has bounded treewidth. These results appear in Section 2. Al- 
though the local-to-global property can be used in conjunction with the algorithm of Elberfeld et al. [10] to 
obtain an improved approximation guarantee for the general setting, we proceed by developing an improved 
0(l/(n|P|)^/^)-approximation algorithm for this problem. Our algorithm is based on a greedy approach that 
employs the local-to-global property in a novel way. The specifics of this algorithm are presented in Section 3. 
We also study two well-motivated variants of the orientation problem, and most notably, show hardness results 
for them. Further details are provided in Section 4. 

2 From Local to Global Orientations 

In this section, we identify a useful structural property of requests crossing through a junction vertex. Infor- 
mally, this property guarantees that if there is an orientation of the local neighborhood of a vertex that locally 
satisfies a set of requests then it can be extended to a global orientation of the complete graph which satisfies 
the same set of requests. Finding a local orientation that maximizes the local satisfiability is a relatively easy 
task, namely, it admits a constant factor approximation algorithm. As a consequence, we can slightly modify 
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the algorithms developed by Elbeifeld et al. [10] so they utilize this property, and obtain improved approxi- 
mation ratios. For example, we eliminate a logarithmic factor from their polylogarithmic approximation ratio 
for the special case that the underlying graph has bounded treewidth. 

We associate each request (s, t) G P with the shortest path p between s and t in the underlying graph. 
Note that in case there are several shortest paths for a request, we associate it with one of them arbitrarily. 
We now introduce some notation and terminology. To better understand the suggested notation, we refer the 
reader to the concrete example in Figure 1 . 

• The local neighborhood of a vertex v is the subgraph that consists of v, all edges incident on v, and 
all vertices adjacent to v. Notice that the local neighborhood graph is a star. 

• Let Pi, be the set of shortest paths of requests that cross v, and let be the corresponding set of local 
paths, that is, the paths of P^, confined to the local neighborhood of v. More precisely, each (global) path 
p ^ Pv gives rise to a (local) path p' G P'^, defined as the intersection of p with the local neighborhood 
of V. Furthermore, for each p' P^, we define its local endpoints s' and t' to be the closest vertices to 
s and t onp that also appear on p', respectively. 

• The local graph orientation problem corresponding to vertex v is defined with respect to the local 
neighborhood graph and the set of local paths P^. The goal is to orient the undirected edges of 

in a way that maximizes the number of satisfied paths in P^. A path is said to be satisfied if there is a 
directed path between its source and target vertices under the orientation. 



Figure 1: (a) Suppose P = {(si, ti), (s2, ^2), (^s, f )} is the set of requests, and note that the shortest paths of 
these requests are marked with the heavy lines. Notice that all these paths cross v. (b) The local neighborhood 
of V, and the corresponding set of local paths. For example, notice that the local endpoints of the request 
(si, ti) are s'^ = vi and t'^ = V2. 

Lemma 2.1. Given an orientation ofGy that satisfies a set of local paths S' C P^ then there is an orientation 
ofG that satisfies the corresponding set of global paths S C P^. 

Proof. We argue that if two local paths p[ , p'2 S S' then the con^esponding global paths pi,p2 S S cannot 
be in conflict. The paths pi and p2 are said to be in conflict if they have a mutual undirected edge that gets 
a different direction when the edges of pi are consistently oriented from its source vertex to its tai^get vertex 
and when the edges of p2 are consistently oriented from its source vertex to its target vertex. Notice that 
establishing this argument completes that proof of the lemma since none of the paths of S can be in conflict 
with another path in S, and therefore, all the paths in S can be simultaneously satisfied by consistently 
orienting each one of them from its source vertex to its tai^get vertex. Note that after one orients those paths, 
the remaining undirected edges of the graph can be oriented in some arbitrary way. 




(a) 



(b) 



3 



For the puipose of establishing the above argument, let us suppose that pi and p2 are in conflict, and attain 
a contradiction. Since pi and p2 are in conflict then there is an undirected edge e = {vi,V2) € Ey that gets 
a different direction when consistently orienting each one of pi and p2 from its source vertex to its target 
vertex. Let us assume without loss of generality that edge e is the closest to v from all conflicting edges. We 
next present a case analysis that depends whether the edge e appears before or after the position of vertex v 
on each of paths pi and p2- Essentially, there are two main cases. To better understand the used notation, we 
refer the reader to the concrete examples in Figure 2. 

Case I: edge e appears after vertex v in both pi and P2- Let us assume without loss of generality that vi 
is closer to v than V2 on pi, and V2 is closer to v than vi on p2. Let di be the distance between v and vi on 
pi, and d2 be the distance between v and V2 on p2. Since pi is a shortest path between si and ti, it must also 
be a shortest path between v and V2. Thus, di + 1 < ^2- Similarly, since p2 is a shortest path between S2 and 
t2, it must also be a shortest path between v and vi, and hence, d2 + I < di. Summing together the above 
inequalities results in di + ^2 + 2 < + (i2, a contradiction. 

We note that the case that the edge e appears before vertex v in both pi and p2 can be handled along the 
same lines with an adjustment to the relative position of v, e.g., the distances need to be defined from vi and 
V2 towards the junction vertex v. 

Case II: edge e appears after vertex v in pi and before vertex v in p2- Let us assume without loss of 
generality that vi is closer to v than V2 on both paths pi and p2. Since p'i,p'2 S S' we know that the edge 
on which pi leaves v and the edge on which p2 enters v must be different. This implies that the subpath 
between v and vi on pi and the subpath between vi and v on p2 are different. Consequently, merging these 
two subpaths creates a cycle in the graph. This contradicts the fact that the graph is acyclic. 

Note that the case that the edge e appears after vertex v in p2 and before vertex v in pi is essentially 
identical to the above case up to a renaming of the paths. ■ 



Figure 2: (a) The case that e appears after v in both pi and p2. (b) The case that e appears after v in pi and 
before v inp2. 

We now concentrate on the computational complexity of the local graph orientation problem correspond- 
ing to a vertex v. One can easily validate that this problem is equivalent to the maximum undirected graph 
orientation problem on a star. Medvedovsky et al. [22] demonstrated that this problem is equivalent to the 
maximum directed cut problem. This latter problem admits constant factor approximation algorithms (see, 
e.g., [11, 20]). In fact, one can easily verify that a random orientation of the undirected edges in the local 
neighborhood satisfies at least 1/4 of the paths of P'^ in expectation. This follows since the maximal length 
of any path in the local neighborhood is at most 2. Furthermore, one can use the method of conditional ex- 
pectations to obtain a detemiinistic orientation that satisfies at least 1/4 of the paths, and consequently, this 
approach is a 1/ 4-approximation for this problem. Combining this result with the local-to-global orientation 
property exhibited in Lemma 2.1 implies the following theorem. 




(a) 



(b) 
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Theorem 2.2. Given a vertex v and a set of requests whose shortest paths cross v, there is a polynomial- 
time algorithm that computes an orientation that satisfies Vl{\P^\) requests. 

We can now modify the algorithms developed by Elbeifeld et al. [10] in accordance with Theorem 2.2, 
and obtain the following improved approximation ratios. We emphasize that the algorithms and their analysis 
follow (up to our modification step) those presented by Elberfeld et al. [10], and thus, we defer them to the 
appendix. The first two theorems present algorithms whose approximation guarantees depend on the treewidth 
and feedback vertex number of the underlying graph. 

Theorem 2.3. There is a polynomial-time algorithm that finds an orientation satisfying Vt{\P\/ {k\ogn)) 
requests when the undirected version of the underlying graph has bounded treewidth k. 

Theorem 2.4. There is a polynomial-time algorithm that finds an orientation satisfying ^{\P\/[k + logn)) 
requests, where k is the minimum number of vertices whose deletion turns the undirected version of the 
underlying graph into a tree. 

We can also improve the approximation ratios of the algorithms presented by Elberfeld et al. [10] for the 
general case, in which there are no structural assumptions on the graph, by a logarithmic factor. 

Theorem 2.5. There is a polynomial-time algorithm that approximates the maximum mixed graph orientation 
problem to within a factor o/ri(l/yA|P|), where A is the maximum length of a shortest source-target path 
in the graph. 

Theorem 2.6. There is a polynomial-time algorithm that approximates the maximum mixed graph orientation 
problem to within a factor ofQ.{l/M^/^), where M = max{n, \P\}. 

Note that we do not provide a proof for the latter theorem since it can be established along the same lines 
of [10], but more importantly, since we next present an algorithm with a better approximation guarantee. 

3 Improved Approximation for The General Case 

In this section, we develop a relatively simple r2(l/(77,|P|)^/'^)-approximation algorithm for the maximum 
mixed graph orientation problem. Our algorithm is based on a greedy approach that employs the local-to- 
global orientation property developed in Section 2. 

The algorithm, formally described below, begins by associating each request (sj, tj) G P with a shortest 
path Pi between Sj and ti in the graph. Then, it greedily orients shortest paths one after the other until all 
the remaining paths are in conflict with many other paths. When this happens, the algorithm concentrates on 
the vertex that is crossed by a maximal number of paths, and utilizes the local-to-global orientation algorithm 
from Theorem 2.2 to complete the orientation of the graph. Recall that two paths pi and p2 are said to be in 
confiict if they have a mutual undirected edge that gets a different direction when the edges of pi and p2 ai"e 
consistently oriented from their source vertex to their target vertex. 

One can easily verify that the algorithm computes a feasible orientation, namely, it assigns a single di- 
rection to each undirected edge. This follows since no conflicting paths are oriented during the main loop 
of the algorithm, and since the algorithm from Theorem 2.2 is known to compute a feasible orientation. We 
next prove that the algorithm satisfies rj(l/(n|P|)^/'^)-fraction of all requests. Clearly, this implies that the 
algorithm achieves (at least) the same approximation guarantee. 

Theorem 3.1. The greedy orientation algorithm satisfies n{l / {n\P\Y/^) -fraction of all requests. 
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Algorithm 1 Greedy Orientation 

Input: A mixed graph G and a collection P ^ZV xV of requests 
Output: An orientation G of G 

1: Let Pi be a shortest path for request {si,ti) G P in G, and let V = U{k} 
2: while there is pi £ V that is in conflict with less than paths in V do 

3: Let Q C "P be the set of paths in conflict with pi 

4: G ^ the graph that results by orienting the edges of pi from Si towards ti in G 

5: V ^V\{QU{p,}) 

6: end while 

7: Let be a vertex that a maximal number of paths in V cross, and let "P^ C "P be that set of paths 
8: G the graph that results by executing the algorithm from Theorem 2.2 with respect to v and P, 

9: return G 



Proof. Let P = |J{P*} '^^e initial collection of shortest paths, and note that |P| = |P|. In addition, let 
P2 C P be the set of paths the remain after the termination of the main loop of the algorithm, and Pi = P\P2. 
Finally, let Ai be the set of paths that our algorithm satisfies during the main loop of the algorithm, and let 
A2 be the set of paths that the algorithm satisfies during the execution of the algorithm from Theorem 2.2. In 
what follows, we prove that \Ai\ = 0(l/(n|P|)i/3) • \Vi\, and \A2\ = $7(l/(n|P|)^/2) • IP2I. Consequently, 
we obtain that the number of paths satisfied by our algorithm is 

The fact that \Ai\ = 0(l/(n|P|)^/^) • \Vi\ easily follows by observing that in each step of the main 
loop of the algorithm, one path is satisfied while less than (n|P|)^/'^ paths are discarded. Hence, we are left 
to prove that \ A2\ = 0(l/(n|P|)^/^) • IP2I. We establish a somewhat stronger result by demonstrating that 
1^2! = f^(l/(n|p2|)i/3) . For this purpose, consider two paths pi,p2 G P2 that are in conflict. We 
associate the conflict between these paths to an arbitrary undirected edge that gets a different direction when 
pi and p2 are oriented, and place one token on this edge. Notice that each path of P2 is in conflict with at 
least (n|P|)^/'^ other paths in P2; otherwise, the main loop would not have terminated. This implies that if 
we place a token for each pair of conflicting paths in P2 as shown before then the undirected edges of G have 
at least (n|P|)V3 . \V2\/2 > n^/''^\V2\*^^ /2 tokens placed on them. As a consequence, there must be a vertex 
that has at least t = IP2 (2n^/^) tokens placed on the undirected edges in its local neighborhood. We next 
ai^gue that if some vertex has t tokens in its local neighborhood then there must be paths that cross that 

vertex. As a result, we attain that the number of paths that cross the vertex v, i.e., the vertex that a maximal 
number of paths from P2 cross, is at least = rj(|P2p/^/n^/'^). By theorem 2.2, our algorithm satisfies 

a constant fraction of these requests, namely, \ A2\ = i}{l/{n\V2\)^^^) ■ |P2|> as required. 

For the purpose of establishing the above argument, consider some vertex u that has t tokens in its local 
neighborhood. Let us focus on some edge e in this local neighborhood that has r paths that traverse in one 
direction and i paths that traverse in the other direction. Notice that such an edge is assigned r • i tokens. This 
implies that if the local neighborhood of u consists only of the edge e then the minimal number of paths that 
cross u corresponds to the solution of min{r + £ : r ■ i = t}. One can easily verify that the solution for this 
expression is r = ^ = \/t, that is, the number of paths is Q{\/t). Note that when there is more than one edge 
in the local neighborhood of u then any path may cross at most two edges. As a result, if we denote the set 
of edges in the local neighborhood of u by Eu, then the minimal number of paths that cross u dominates the 
solution of min{^gg^^(re + i?e)/2 : Ylei^e'^e) = t}; here, rgand^e indicate the number of paths traversing 
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edge e in one direction and the other direction, respectively. One can easily demonstrate that the solution for 
the above expression is obtained by assigning non-zero values only to one pair of , £e variables, namely, it 
is equivalent to the solution for the single edge case. ■ 



4 Other Orientation Variants 

In this section, we study two well-motivated variants of the orientation problem: the first is maximum mixed 
graph orientation with fixed paths, and the other is maximum mixed grid orientation. 

4.1 Orientation with fixed paths 

We consider the maximum mixed graph orientation with fixed paths problem. This variant is identical to the 
maximum mixed graph orientation problem with the exception that each request (s, t) G P is also associated 
with a fixed path p from s to t in the graph. With this modified definition in mind, a request {s, t) is satisfied 
only if the edges of the path p are oriented from the vertex s towards the vertex t. Note that this variant 
is seemingly simpler than maximum mixed graph orientation since the only computational task is to decide 
which requests to satisfy, and there is no need to decide which paths will be used to satisfy those requests. 
This is also one of our motivations for studying this variant, hoping that it will shed some light on the original 
problem that would lead to a reduction in the gap between its lower and upper approximation bounds. 

We prove that the maximum mixed graph orientation with fixed paths problem is NP-hard to approximate 
to within a factor of maxll/lPl^"*^, l/m^/^"*^}, for any e > 0. In fact, we establish this result even when 
the underlying graph is undirected. As a consequence, we attain that this problem is provably harder than the 
maximum mixed (or undirected) graph orientation problem, although it may seem simpler at first glance. Our 
proof is based on showing that the problem under consideration captures the well-known maximum indepen- 
dent set problem as a special case. 

A hardness of approximation result. An input instance for the maximum independent set problem consists 
of an undirected graph G' = {V',E'). The goal is to find an independent set of maximum size in the graph. An 
independent set is a collection of vertices that do not have any edges between them. This problem is known to 
be NP-hard to approximate within a factor of l/jy'l^"*^, for any e > [25]. We next show a value-preserving 
reduction from this problem to our maximum mixed graph orientation with fixed paths problem. 

Given an input instance of maximum independent set, we construct an input instance for our problem that 
consists of the undirected graph presented in Figure 3(a). Specifically, we begin by creating a graph with 
n' = \V'\ pairs of Si,ti vertices, corresponding to the vertices of G', such that each such pair is connected 
by a path pi. We intersect all these connecting paths in a grid-like fashion. Then, each intersection point is 
replaced by one of the gadgets exhibited in Figures 3(b) and 3(c). The gadget g-ij that replaces the intersection 
point of paths pi and pj has 4 vertices: Vi and Ui that are appropriately added to path pi, and vj and Uj that 
are appropriately added to path pj. The edges within the gadget has the form described in Figure 3(b) if 
(z, j) ^ E', or the form described in Figure 3(c) if G E'. In the latter case, the path pj is also modified 
to consist of the vertices Vi and Ui, so its subpath inside the gadget is {vj,Ui,Vi,Uj). In addition, the set of 
requests P for our problem consist of all n' pairs (sj, tj) with their corresponding path pi. 

One can easily validate that a solution S C V' for the maximum independent set problem implies an 
orientation in the newly-created instance that satisfies the same number of requests. Specifically, if i,j G S 
then clearly (i, j) ^ E', and thus, the paths pi and pj do not share edges. As a result, one can simultaneously 
satisfy both request {si,ti) and {sj,tj) by orienting each of their paths from its source vertex to its target 
vertex. Conversely, it is not difficult to verify that given an orientation in newly-created instance that satisfies 
some set of requests, one can perform a similar value-preserving transformation in the opposite direction. In 
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ti t2 tn' Uj 



(a) (c) 

Figure 3: (a) The graph resulting from the reduction, (b) The gadget that is used in case {i, j) ^ E'. Note that 
the edges {vi, ui) and {vj,Uj) do not intersect, (c) The gadget that is used in case E E'. 

particular, notice that if i ^ S then there must be j G 5 such that (i, j) e E', and hence, the paths pi and 
Pj require to orient the edge {vi,u,i) in the gadget gij in conflicting directions. Consequently, one cannot 
simultaneously satisfy both underlying requests. 

As a result of this value-preserving reduction, and in conjunction with the hardness result presented by 
Zuckerman [25], we attain the following inapproximability result. Recall that |P| = \ V'\, and notice that the 
number of edges in the newly-created instance is m = 0(|F'p). 

Theorem 4.1. The maximum mixed graph orientation with fixed paths problem is NP-hard to approximate 
within afactorofmax{l/\P\^~'^, l/m^^'^"''}, for any e > 0. 

4.2 Orientation in grid networks 

We study the maximum mixed grid orientation problem. This variant is identical to the maximum mixed 
graph orientation problem with the additional restriction that the graph is a grid. A n x m grid network is a 
graph with a vertex set F = {1, . . . , n} x {1, . . . , m}, and an edge set E consisting of horizontal edges, i.e., 
edges + 1)) for all j = {1, . . . ,m — 1}, and vertical edges, i.e., edges {i + for all 

i = {1, . . . , n — 1}. Note that the study of this variant is motivated by applications in networking. 

We prove that the maximum mixed grid orientation problem is at least as hard as the maximum directed 
cut problem. Consequently, approximating our problem within factors of 12/13 « 0.923 and acw ~ 0.878 
is NP-hard and Unique Game-hard, respectively. Interestingly, this finding comes in contrast with the results 
attainable for the undirected grid setting. This latter setting can be solved to optimality in polynomial-time, 
and in particular, when the grid is not a path, that is, when n,m > 1, all the requests in P can be satisfied. 

A hardness of approximation result. An input instance for the maximum directed cut problem consists of a 
directed graph G' = {V', E'). The goal is to find a directed cut of maximum size in the graph. The size of 
a cut A C y is the number of directed edges {u, v) G E' such that u ^ A and v £ V' \ A. Approximating 
this problem within factors of 12/13 « 0.923 and acw ~ 0.878 is known to be NP-hard [18] and Unique 
Games-hard [19], respectively. In what follows, we present a value-preserving reduction from this problem to 
our maximum mixed grid orientation problem. 
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Given an input instance of maximum directed cut, we construct an input instance for our problem which 
consists of the mixed grid presented in Figure 4(a). Specifically, we create a grid whose dimensions are 
n = 2\V'\ — 1 and m = 3. We associate each vertex Vi € V with the vertex (2i — 1, 1) in the grid. The 
edges incident on each vertex {2i, 1) in the grid are oriented away from that vertex, and the edges along the 
perimeter of the sub-grid that consists of the second and third vertex columns are oriented in a way that creates 
a directed cycle. In addition, the set of requests for our problem is defined to be P = E'. 

One can validate that a solution A V for the maximum directed cut problem implies an orientation 
in the newly-created instance that satisfies the same number of requests. Specifically, if Vi ^ A we orient 
the single undirected edge incident on vertex Vi of the grid away from that vertex, and if G V \ A we 
orient that edge towards vertex Vi. Then, it is easy to see that if an edge of E' is cut by the solution A 
then the corresponding request is satisfied in the orientation. Conversely, it is not difficult to verify that 
given an orientation in newly-created instance that satisfies some set of requests, one can perform a similar 
value-preserving transformation in the opposite direction. In particular, this side of the proof builds upon the 
observation that any request {vi, Vj) may only be satisfied by a path that crosses the undirected edges incident 
on Vi and Vj. The orientation of all those undirected edges define the cut in the initial problem. 

As a result of this value-preserving reduction, and in conjunction with the hardness results presented by 
Hastad [18] and Khot et al. [19], we attain the following inapproximability result. 

Theorem 4.2. The maximum mixed grid orientation problem is NP-hard to approximate within a factor of 
12/13 ~ 0.923, and Unique Games-hard to approximate within a factor of uq^j^ ~ 0.878. 
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Figure 4: (a) The grid resulting from the reduction. Note that n' = \V'\. (b) An orientation of an undirected 
grid that admits a directed path between any two vertices. 



Orientation of undirected grids. The above-mentioned hardness result comes in contrast with the results 
attainable for the undirected grid setting. This latter setting can be solved to optimality in polynomial-time. 
Specifically, when the grid is a path, i.e., when either m or n equals 1, there are optimal polynomial-time 
algorithms for the problem [22, 8], and when n,m > 1, there is a simple orientation that satisfies all the 
requests in P. This orientation can be obtained by creating a directed cycle along the perimeter of the grid, 
and then, orienting all the remaining horizontal and vertical edges consistently. A concrete example of such 
an orientation is presented in Figure 4(b). One can easily prove that this orientation admits a directed path 
between any two vertices of the graph. 



9 



References 

[1] Y. Afek and A. Bremler-Barr. Self-stabilizing unidirectional network algorithms by power supply. 
Chicago J. Theor. Comput. Sci., 1998. 

[2] Y. Afek and E. Gafni. Distributed algorithms for unidirectional networks. SIAMJ. Comput., 23(6):1152- 
1178, 1994. 

[3] E. M. Arkin and R. Hassin. A note on orientations of mixed graphs. Discrete Applied Mathematics, 
116(3):27 1-278, 2002. 

[4] V. Bafna, P. Berman, and T. Fujito. A 2-approximation algorithm for the undirected feedback vertex set 
problem. SIAMJ. Discrete Math., 12(3):289-297, 1999. 

[5] A. Becker and D. Geiger. Approximation algorithms for the loop cutset problem. In Proceedings 10th 
Annual Conference on Uncertainty in Artificial Intelligence, pages 60-68, 1994. 

[6] H. L. Bodlaender. A linear-time algorithm for finding tree-decompositions of small treewidth. SIAM J. 
Com/JMf., 25(6):1305-1317, 1996. 

[7] F. A. Chudak, M. X. Goemans, D. S. Hochbaum, and D. P. Williamson. A primal-dual interpretation of 
two 2-approximation algorithms for the feedback vertex set problem in undirected graphs. Oper Res. 
Lf-ff., 22(4-5): 11 1-1 18, 1998. 

[8] B. Dorn, F. Hiiffner, D. Kriiger, R. Niedermeier, and J. Uhlmann. Exploiting bounded signal flow for 
graph orientation based on cause-effect pairs. In Proceedings 1st International ICST Conference on 
Theory and Practice of Algorithms in (Computer) Systems, pages 104-115, 2011. 

[9] M. Elberfeld, V. Bafna, I. Gamzu, A. Medvedovsky, D. Segev, D. Silverbush, U. Zwick, and R. Sharan. 
On the approximability of reachability-preserving network orientations. Internet Mathematics, 7:209- 
232, 2011. 

[10] M. Elberfeld, D. Segev, C. R. Davidson, D. Silverbush, and R. Sharan. Approximation algorithms for 
orienting mixed graphs. In Proceedings 22nd Annual Symposium on Combinatorial Pattern Matching, 
pages 416^28, 2011. 

[11] U. Feige and M. X. Goemans. Aproximating the value of two prover proof systems, with applications 
to MAX 2SAT and MAX DICUT. In Proceedings 3rd Israel Symposium on Theory and Computing 
Systems, pages 182-189, 1995. 

[12] S. Fields. High-throughput two-hybrid analysis: The promise and the peril. The FEBS Journal, 
272(21):5391-5399, 2005. 

[13] G. N. Frederickson and D. B. Johnson. Generating and searching sets induced by networks. In Pro- 
ceedings 7th International Colloquium on Automata, Languages and Programming, pages 221-233, 
1980. 

[14] I. Gamzu and D. Segev. A sublogaiithmic approximation for highway and toUbooth pricing. In Pro- 
ceedings 37th International Colloquium on Automata, Languages and Programming, pages 582-593, 
2010. 



10 



[15] I. Gamzu, D. Segev, and R. Shaian. Improved orientations of physical networks. In Proceedings 10th 
International Workshop on Algorithms in Bioinformatics, pages 215-225, 2010. 

[16] A. C. Gavin, M. Bosche, R. Krause, R Grandi, M. Marzioch, A. Bauer, J. Schultz, J. M. Rick, A. M. 
Michon, C. M. Crucial, M. Remor, C. Hofert, M. Schelder, M. Brajenovic, H. Ruffner, A. Merino, 
K. Klein, M. Hudak, D. Dickson, T. Rudi, V. Gnau, A. Bauch, S. Bastuck, B. Huhse, C. Leutwein, M. A. 
Heurtier, R. R. Copley, A. Edelmann, E. Querfurth, V. Rybin, G. Drewes, M. Raida, T. Bouwmeester, 
R Bork, B. Seraphin, B. Kuster, G. Neubauer, and G. Superti-Furga. Functional organization of the yeast 
proteome by systematic analysis of protein complexes. Nature, 415:141-147, 2002. 

[17] S. L. Hakimi, E. F. Schmeichel, and N. E. Young. Orienting graphs to optimize reachability. Information 
Processing Letters, 63(5):229-235, 1997. 

[18] J. Hastad. Some optimal inapproximability results. Journal of the ACM, A%{A):19%-%59 , 2001. 

[19] S. Khot, G. Kindler, E. Mossel, and R. O'Donnell. Optimal inapproximability results for MAX-CUT 
and other 2-variable CSPs? SIAM Journal on Computing, 37(l):319-357, 2007. 

[20] M. Lewin, D. Livnat, and U. Zwick. Improved rounding techniques for the MAX 2-SAT and MAX 
DI-CUT problems. In Proceedings 9th International Conference on Integer Programming and Combi- 
natorial Optimization, pages 67-82, 2002. 

[21] M. K. Marina and S. R. Das. Routing performance in the presence of unidirectional links in multihop 
wireless networks. In Proceedings 3rd ACM Interational Symposium on Mobile Ad Hoc Networking and 
Computing, pages 12-23, 2002. 

[22] A. Medvedovsky, V. Bafna, U. Zwick, and R. Sharan. An algorithm for orienting graphs based on 
cause-effect pairs and its applications to orienting protein networks. In Proceedings 8th International 
Workshop on Algorithms in Bioinformatics, pages 222-232, 2008. 

[23] D. Silverbush, M. Elberfeld, and R. Sharan. Optimally orienting physical networks. In Proceedings 
15th Annual International Conference on Research in Computational Molecular Biology, pages 424- 
436, 2011. 

[24] C.-H. Yeang, T. Ideker, and T. Jaakkola. Physical network models. Journal of Computational Biology, 
ll(2/3):243-262, 2004. 

[25] D. Zuckerman. Linear- degree extractors and the inapproximability of max clique and chromatic number. 
Theory of Computing, 3:103-128, 2007. 

A Additional Details 

In this section, we complete the details omitted from the main part of the paper. 
A.l Proof of Theorem 2.3 

The algorithm begins by computing a tree decomposition of width k for the undirected version of the under- 
lying graph. Note that this task can be done in polynomial-time since our graph has bounded treewidth [6]. 
A tree decomposition is a pair (T, X), where X = {Xi, . . . , X^} is a collection of subsets such that each 
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Xi C V, and T is a tree whose nodes are the subsets in X. Note that i = poly(n) in our case. The decompo- 
sition satisfies the following properties: (1) [jXi = V, (2) the incident vertices of every edge of the graph are 
contained in some subset in X, and (3) if Xi and Xj contain a vertex v then all the nodes X^ in the unique 
path between Xi and Xj contain v as well. The width of the tree is defined to be max |Xj| — 1. 

Given the tree decomposition (T, X), the algorithm proceeds by computing a shortest path for each re- 
quest in P. The paths are then classified into at most [log^] = O(logn) classes such that for every class, an 
orientation that satisfies r2(l/A;)-fraction of its paths can be efficiently computed. As a consequence, by sep- 
ai^ately computing an orientation for each class, and then picking the option that satisfies the highest number 
of paths, we ai^e guaranteed to satisfy Q{\P\/{k log n)) of all the requests. 

For the purpose of constructing the first class, we find a centroid node Xt of T, that is, a node whose 
removal breaks the tree into a collection of subtrees, each of which has at most half of the vertices in T. Note 
that any tree has a centroid (see, e.g., [13]). We assign all the paths that cross a vertex from Xt = {vi, . . . , Vr} 
to class Ci. We further partition Ci into r collections Ci^i, . . . ,Ci^r such that a path p is assigned to the 
collection Cij if it crosses vj but does not cross any of the vertices in {vi, . . . , vj^i}. Notice that we can 
satisfy r2(|Cij|) paths from the collection j by applying Theorem 2.2. One can now easily validate that 
executing the mentioned algorithm on each collection separately, and then picking the option that satisfies the 
highest number of paths results in an orientation satisfying ^{\Ci\/k) requests since r < A; + 1. 

To construct the second class, we first remove the node Xt from T to obtain a forest of tree decomposi- 
tions. For each tree decomposition, we compute a centroid node, and in the same way as above, we assign 
a path to C2 if it crosses a vertex from the subsets associated with these centroid nodes. Note that we only 
assign paths that were not assigned to the first class. Using the same arguments as above, we can compute 
an orientation that satisfies 0(|C2|/fc) requests. In particular, one can validate that each path crosses vertices 
from exactly one centroid node; otherwise, it should have been assigned to the first class by properties (2) and 
(3) of the tree decomposition. We now proceed recursively in the same way to construct the other classes as 
long as the decompositions under consideration are not empty. Since the maximal size of a subtree decreases 
by at least half in each level of the recursion, this process terminates within [log f\ steps, and hence, there ai^e 
indeed at most [log i] classes. 

A.2 Proof of Theorem 2.4 

The algorithm begins by finding a feedback vertex set F = {vi, . . . ,Vi} in the undirected version of the 
graph, namely, a set of vertices whose removal turns the underlying undirected graph into a tree. Although 
the computational task of finding a feedback vertex set with a minimum cardinality is NP-hard, there is a 
2-approximation algorithm for this problem [5, 4, 7]. Therefore, we may assume that the cardinality of that 
set satisfies i < 2k. The algorithm proceeds by computing a shortest path pi for each request {si,ti) G P. 
Then, each path pi is classified into one of ^ + 1 classes: if pi crosses the vertex Vj and none of the vertices in 
{vi, . . . , Vj^i} then it is assigned to class Cj; otherwise, if pi does not cross any of the vertices of F, then it 
is assigned to class C^+i. Notice that we can satisfy r2(|Cj |) paths from any class j by applying Theorem 2.2. 
Also notice that by deleting the vertices of F from the graph G, we obtain a mixed graph which is a forest 
of trees, and all the paths in Cc^i still remain connected. This mixed tree orientation setting is known to 
admit an efficient ^^(1/ log n) -approximation algorithm [10], and thus, we can satisfy f^dC^+i |/ log n) paths 
from the class i + 1. One can now easily validate that executing the mentioned algorithms on each class 
separately, and then picking the option that satisfies the highest number of paths results in an orientation 
satisfying i}{\P\/{k + log??,)) of all requests. 
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A.3 Proof of Theorem 2.5 



The algorithm computes a shortest path pi for each request (sj,ti) € P. Then, it considers those shortest 
paths in some arbitrary order, and orients them one after the other. In paiticulai^, when a path p is oriented, 
all the pending paths that are in conflict with p aie discai^ded. This greedy orientation procedure continues 
as long as a path p under consideration is not in conflict with more than yA|P| pending paths. When this 
happens, there must be some vertex v on p that at least ^A|P|/A = -y/lPl/A pending paths cross. This 
claim holds since the length of p is known to be at most A. The algorithm then employs the local-to-global 
orientation algorithm from Theorem 2.2 with respect to the vertex v and the corresponding set of pending 
paths to complete the orientation of the graph. 

One can easily verify that the algorithm computes a feasible orientation. Therefore, we next prove that the 
algorithm satisfies A|P|)-fraction of all requests. Notice that if the algorithm employs our local-to- 

global orientation algorithm then the resulting orientation indeed attains the desired performance guarantee 
since the number of paths satisfied only by this step is Q{^y\P\/ A) = ^}{l/^^/A\P\) ■ \P\. Hence, we 
may assume that the algorithm only makes greedy orientation steps. In this case, the resulting orientation 
clearly achieves the desired performance guarantee since one path is satisfied while at most yA|P| paths are 
discarded in each of those greedy steps. 
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